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The NASA initiative to collect an asteroid, the Asteroid Robotic Redirect Mission (ARRM), 
is currently investigating the option of retrieving a boulder from an asteroid, demonstrat- 
ing planetary defense with an enhanced gravity tractor technique, and returning it to a lunar 
orbit. Techniques for accomplishing this are being investigated by the Satellite Servicing 
Capabilities Office (SSCO) at NASA GSFC in collaboration with JPL, NASA JSC, LaRC, 
and Draper Laboratory, Inc. Two critical phases of the mission are the descent to the boulder 
and the Enhanced Gravity Tractor demonstration. A linear covariance analysis is done for 
these phases to assess the feasibility of these concepts with the proposed design of the sensor 
and actuator suite of the Asteroid Redirect Vehicle (ARV). The sensor suite for this analysis 
includes a wide field of view camera, LiDAR, and an IMU. The proposed asteroid of interest 
is currently the C-type asteroid 2008 EV5, a carbonaceous chondrite that is of high interest 
to the scientific community. This paper presents an overview of the linear covariance anal- 
ysis techniques and simulation tool, provides sensor and actuator models, and addresses the 
feasibility of descending to the surface of the asteroid within allocated requirements as well 
as the possibility of maintaining a halo orbit to demonstrate the Enhanced Gravity Tractor 
technique. 

INTRODUCTION 

The Asteroid Robotic Redirect Mission is evaluating two different options to fulfill the primary goal of 
redirecting a small asteroid or boulder and inserting it into a lunar orbit. Option A is to target and capture 
a ten-meter class asteroid (up to 1000 metric tons,) whereas Option B is to target a larger asteroid, collect a 
two- to four-meter boulder (up to 70 metric tons) from the surface, and perform a planetary defense demon- 
stration. 1 Techniques for accomplishing this are being investigated by the Satellite Servicing Capabilities 
Office (SSCO) at NASA GSFC in collaboration with JPL, NASA JSC, LaRC, and Draper Laboratory, Inc. 
The planetary defense phase consists of an Enhanced Gravity Tractor (EGT) demonstration, which uses the 
gravitational force of the spacecraft and boulder to deflect the orbit of the asteroid. This concept was first 
proposed by Lu and Love in 2005. 2 A follow-on mission is being planned to then send astronauts to study the 
asteroid while it is in orbit around the moon. Previous work was done to develop a similar mission concept for 
Itokawa. 3 The focus of this paper will be to assess the feasibility of Option B given the current configuration 
for asteroid 2008 EV5. 

Two critical phases of the mission are the descent to the surface of the asteroid and the planetary defense 
demonstration. A linear covariance analysis is performed for each of these phases. The next section gives a 
brief overview of linear covariance (LinCov) analysis and the tool utilized in this research. Then an overview 
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of the dynamics as well as the sensor and actuator models is presented. Finally the concept of operations, 
control logic, and results are discussed for each flight phase. 

LINEAR COVARIANCE ANALYSIS OVERVIEW 

In general, to perform the necessary analysis and evaluation of a GN&C system, there are several key 
variables of interest, including environment dispersions Sx, navigation dispersions Sx., actual navigation 
error Se, and onboard navigation error Se. These variables characterize the system’s performance and are 
used to develop and validate mission objectives and requirements. 4 Extensive effort and resources are often 
allocated for the very purpose of producing and analyzing the quantities of these system parameters. For this 
paper, they are instrumental in characterizing the overall system and navigation performance for descending 
to the asteroid’s surface and demonstrating the enhanced gravity tractor concept. This section is dedicated 
to formally defining these metrics and summarizing the analysis techniques and simulation tools adopted for 
obtaining them. 

Performance Metric Definitions 

The environment dispersions are defined as the difference between the environment state x and the nominal 
state x. The environment state is an (n x 1) vector that represents the true or actual state. The nominal state 
is also an (n x 1) vector that represents the desired or reference state. The covariance of the environment 
dispersions, D, indicates how precisely the system can follow a desired trajectory. 

Sx = x — x, D = E [£x£x t ] (1) 

The environment (or true) dispersions are often referred to as simply dispersions (e.g., trajectory dispersions, 
position dispersions, relative dispersions). 

The navigation dispersions are defined as the difference between the navigation state x and the nominal 
state. The navigation state is an (n x 1) vector that represents the estimated state. 

Sx = it — Nx, D = E [£x£x t ] (2) 

The matrix N is an (h x n) mapping matrix that defines the estimated state in terms of the true and nominal 
state. It typically cancels the attitude rate state when gyro measurements are incorporated in lieu of an angular 
rate estimate. An inverse mapping takes the estimated state to the true state, represented with an (n x n) matrix 
N T . The covariance of the navigation dispersions, D, reflect how precisely the onboard system thinks it can 
follow a prescribed reference trajectory. 

The true navigation error is the difference between the environment and navigation states. It is also the 
difference between the environment and the navigation dispersions. 

Se = Nx — x = N Sx — Sx , P — E [£e£e T ] (3) 

The covariance of the true navigation error, P, characterizes how precisely the onboard navigation system can 
determine the actual state. The onboard navigation error is never computed but is used to develop the onboard 
navigation filter equations. It is defined as the difference between the design state, x, and the navigation state. 

Se = x — x, P = E [£e£e T ] (4) 

The covariance of the onboard navigation error, P, shows how precisely the onboard navigation system 
thinks it can determine the actual state. The performance of the onboard navigation system is determined by 
comparing P to the actual navigation performance P. It is the covariances of the true dispersions, navigation 
dispersions, true navigation error, and the onboard navigation error that are ultimately used to analyze and 
assess the performance of a proposed GN&C system. 
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Obtaining the Performance Metrics 

A common approach to obtain these performance metrics is to use a Monte Carlo simulation as shown in 
Figure 1 , where the sample statistics of hundreds or thousands of runs are used to numerically compute the 
desired covariance matrices. 


D = 


D 


N — 




This same statistical information can be obtained using linear covariance analysis techniques 5,6,7,8 by 
directly propagating, updating, and correcting an augmented state covariance matrix C, 

C = E [£X£X t ] (5) 

where the augmented state SX. T = [£x T £x T ] consists of the true dispersions and the navigation dispersions. 
Notice that by simply pre- and post-multiplying the augmented state covariance matrix by the following 
matrices, the covariance matrices for the trajectory dispersions, navigation dispersions, and the navigation 
error can be obtained. 
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Consequently, Monte Carlo and linear covariance analysis techniques provide a complimentary analysis pack- 
age since they each generate the same statistical information using different approaches. The strengths and 
weaknesses of one technique often offset the other. Such is the case when developing and analyzing the 
preliminary trajectory design and concept of operations for close-proximity asteroid operations. 



Figure 1. GN&C Performance Metrics in a Generic Monte Carlo Simulation 4 


SIMULATION OVERVIEW 

For this preliminary design phase and feasibility study, a linear covariance analysis tool is utilized. Minor 
adjustments are made to a previously developed rendezvous and docking LinCov simulation 8, 9 by including 
the gravitational attraction between both the target (asteroid) and chaser (ARV) bodies, which is typically 
neglected in rendezvous scenarios between two orbiting vehicles. For this simulation, the absolute states of 
the ARV and the target asteroid are modeled in the Sun-Centered Inertial reference frame. The asteroid 2008 
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EV5 is considered in this paper. It has a mass of rat=1.05E+l 1 kg, 10 a gravitational parameter of 7 m 3 /s, and 
a radius of 200 m. The asteroid trajectory and initial position and velocity is obtained from a Horizons .bsp 
file. 11 The initial attitude is q t = [ -0.1812, -0.9587, 0.1929, 0.1039], and the asteroid rotates about its north 
pole with a four hour period. The asteroid inertia, J t , has moments of inertia [4.42855, 4.48630, 4.76455] 
10 15 kg-m 2 computed using dimensions from Reference 10 and zero products of inertia. 

The ARV mass is m c - 9000 kg, and ARV inertia, J c , has moments of inertia [81886.5, 120382.5, 117450] 
kg-m 2 and zero products of inertia. The simulation starts on December 31, 2023, and uses a fourth-order 
Runge-Kutta method with a 1 s time step. A spherical gravity field and J2000 reference frame centered at the 
Sun are utilized.The relative initial state in Asteroid-Centered Inertial reference frame is a position of [-152.85 
-158.54 333.91] m and a velocity of [-7.161 85.66 37.397] mm/s with and initial attitude of q c =[0.1988, 
-0.0032, -0.9578, -0.2076]. The relative states are estimated using a suite of sensors described later in this 
section. The relative 3a uncertainty in the initial states is as follows: 10 m in position, 5 mm/s in velocity, 
1 deg in attitude, 1 arcsec/s in attitude rate. The 3 a process noise is 1.1727E-6 m/s/ a/s for translation and 
0.00 IE-9 rad/ </ s for rotation. The translational and rotational maneuver execution errors are given in Tables 
4 and 5. The concept of operations for each phase will be discussed in more detail later. 

Dynamic Modeling 

The asteroid rendezvous simulation models include the chaser and target translation and rotational dynam- 
ics; disturbance accelerations and torques; strap-down inertial measurement units (IMUs), optical navigation 
camera, and FPOSE (from Flash LiDAR pose) measurements; and the torques and velocity changes gener- 
ated by momentum wheels and thrusters. The truth model state is defined by 13 target states x t , 13 chaser 
states x c , and parameter states associated with various aspects of the sensors used onboard the chaser x p . 
The 13 asteroid target states include the inertial position and velocity vectors, the quaternion defining the 
orientation of the target with respect to the inertial frame, and the target’s angular rate coordinatized in the 
target reference frame. Similarly, the 13 chaser states are inertial position and velocity, the inertial-to-body 
quaternion, and angular rate. The k parameter states such as the sensor or actuator biases, misalignments, 
scale factor terms, and other error sources are modeled as first-order Markov processes with time constants, 
Tk- 

x = [x t ; x c ; x p ] T (6) 


x t = [r t ; v t ; qf; 


v c ; q*; 


[xpl, Xp2) •••! Xpfc] 


(7) 


The dynamics for the chaser, target, and parameter states are defined as, 


ft 

= Vt 

(8) 

Ac 

v c 

Vt 

= S t + T] at 

(9) 


<3* 

= 

(10) 

• c 


= J t _1 [ T S - w t x Jf^t] +v at 

(ID 

(jj 

^ c 

X-pk 


= v c (12) 

= g c + a ctrl + Va c 

= q? (14) 

= 3 c l [ T ctrl + TC g-U C c X J C <<] + ViP) 
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where and J c are the target and chaser inertia matrices, respectively. The accelerations due to gravity 
acting on the target and chaser vehicles, g t and g c , are based on multi-body point-mass gravity terms. The 
gravitational torques, r l g and r c g , are derived from Wertz. 12 The control inputs, a ctr i and r c ctrl , are the 
accelerations and torques executed by the actuators on the chaser spacecraft defined in Eq. (27) and Eq. (29) 
respectively. The random disturbances, 77 , 77 , rj a , and r \ OLc , are included in the models to account for 

unmodeled forces and torques acting on each body such as drag, solar radiation pressure, venting gases, etc. 


Sensor Modeling 

The analysis assumes several different sensors, including an optical camera, FPOSE, and an IMU. The 
mounting configurations of the cameras are illustrated in Figure 2. 
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Figure 2. ARV and Sensor Locations 


Measurement Frequency 

10 

s 

Distortion Coefficient 

[0.06 0 0] 


Resolution 

2592 x 1944 


Focal Length 

4.7 

mm 

Medium Field of View (MFOV) 

10 

deg 

Wide Field of View (WFOV) 

30 

deg 

Measurement Noise 

1 

pixels 

Pixel Density 

2758 

pixels/mm 


Table 1. Camera Parameters 


Optical Camera Model The optical camera model is based on the models described in the references by 
Owen 13 and D’ Souza. 14 The parameters for the cameras currently baselined for the mission and used in 
this simulation are given in Table 1. A process called Stereophotoclinometry (SPC) 15 will be used for the 
optical navigation. This process includes generating maplets and correlating images, and landmarks within 
images, with those maplets. Landmark locations are identified as columns of the matrix l[ mr/c . The landmark 
positions are specified in the target body frame and transformed to the inertial frame as ll mrk = TJl[ mr/c . 
The location of the landmark in the inertial frame is the position of the asteroid plus the landmark position 
relative to the center of the asteroid r l lmrk = r\ + 1 \ mrk - Assuming the camera is at the center-of-mass of the 
chaser vehicle, if the chaser position is known in the inertial frame r l c then the apparent vector in the inertial 
frame is simply the difference a l lmrk = C lrnrk — r l c , as depicted in Figure 3(b). 




(b) Apparent position of land- 
mark, a, position of camera r c , 
position of asteroid r t, and po- 
sition of landmark with respect 
to the asteroid limrk 


Figure 3. Optical Navigation Model Definitions 

The apparent location of the landmark feature relative to the camera in the optical camera frame, a° mr/c , as 
shown in Figure 3(a) is 

A imrk = A. ss [M N L] T°a\ mrk (17) 

where T? is the rotation matrix from the inertial frame to the optical camera frame. Expanding Eq. (17), the 
apparent location of the landmark features can be expressed in terms of the navigation states, which include 
the target inertial position rj, the target orientation uncertainty 0\, the chaser inertial position r£, and the 
chaser attitude error 6 C C . 


x lmrk 


= T°(I - [0“x])T ? v\ + Tj(I + $x])lf mrfc - 4 


( 18 ) 


At the detector, the image appears inverted with pixel coordinates (x,y) which are a function of the apparent 
location of the landmarks A = a° mr/e , the camera focal length /, and the pixel biases, b x and b y . 
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The logic to determine whether the landmark is visible is then executed. Radial distortion is added which dis- 
torts the image toward or away from the optical axis. Tip and tilt terms are added to account for misalignment 
due to the detector not being perfectly perpendicular to the optical axis. Then the corrected image locations 



are transformed to pixel, s and l line coordinates, 



where K contains the reciprocal of the pixel dimensions. 13 Additional logic has been implemented to check 
whether each landmark is within the FOV of the camera and to enforce a minimum maplet resolution of 10 
cm. 

FPOSE Model FPOSE uses LiDAR data to generate a six degree-of-freedom position and orientation 
measurement relative to a target object. This will be used for the descent phase and the target will be the 
boulder of interest. The measurement is modeled with a constant and a Markov bias. 16 The relative position 
vector in the FPOSE frame is given by the equation 

Hel = T b T i ( r * - r c) + tf + ft p + T] f p , (22) 

where r t and r c are the target and chaser inertial position vectors, b f is the relative position constant bias, 
(3* is the relative position Markov bias, and r]f is the relative position measurement noise. 

~f 

The relative attitude measurement in the FPOSE sensor frame is processed as a derived measurement, G rel . 
It is effectively the residual to be processed by the onboard navigation filter, 

I-[e f rel x]= T{ [flflf}] (23) 

where the estimate of the derived relative attitude measurement is a function of the target attitude uncertainty 
Q\, the chaser attitude error state Q c c , the constant bias , a Markov bias f3g, and noise tJq. 

Hel = T {e\- 0 c c + b f e +(3 f e + rjg (24) 

The Markov and constant bias 3cr uncertainties for translation and orientation were 57 cm and 1 deg with 
a time constant of 600 s and 19 cm and 1 degree of measurement noise. The FPOSE uncertainty parameters 
used for this simulation were obtained from test data with a satellite mock-up as the target. A test campaign 
to collect data that more closely represents this scenario is planned for the future. 

IMU Model The gyro model is based upon a package of three orthogonal strapdown gyros, each measur- 
ing the chaser’s angular velocity along its input axis. The measured angular velocity is a function of the true 
angular rate u? c with a gyro bias P % ™ , scale factor s^ mn , misalignment fjd™, and the gyro angular random 
walk i^ mu . 

u imu = T (*4 mw ) {I + [s^ mn \] } T^ u u c + ft™ + (25) 

The gyro performance parameters are summarized in Table 2. 

The measured accelerations from the accelerometers are a function of the non-gravitational accelerations 
a\ the accelerometer bias P % ™, scale factor s l J nu , misalignment fi l ™ , and the accelerometer random walk 

^jimu 

a imu = T ) {I + [sr M \] } Ti mu T? a* + ft™ + vft u (26) 

The accelerometer performance parameters are summarized in Table 3. 
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Constant Bias 


Constant Bias 


Angle Random Walk Uncertainty 

0.018 

deg/\/hr (3cr) 

Bias Uncertainty 

0.45 

deg/hr (3cr) 

Scale Factor Uncertainty 

15 

ppm (3cr ) 

Misalignment Uncertainty 

75 

arcsec (3cr) 


Markov Bias 


Time Constant 

100 

s 

SS Uncertainty 

0.09 

deg/hr (3(j) 


Markov Scale Factor 


Time Constant 

100 

s 

SS Uncertainty 

100 

ppm (3cr ) 


Table 2. Gyro Parameters 


Velocity Random Walk Uncertainty 

150 

lig/Vni (3a) 

Bias Uncertainty 

300 

Mg (3cr) 

Scale Factor Uncertainty 

525 

ppm (3cr ) 

Misalignment Uncertainty 

45 

arcsec (3cr) 


Markov Bias 


Time Constant 

60 

s 

SS Uncertainty 

45 

Mg (3cr) 


Markov Scale Factor 


Time Constant 

100 

s 

SS Uncertainty 

525 

ppm (3a) 


Table 3. Accelerometer Parameters 


Actuators 


Actuators are used to manipulate the spacecraft’s translational and rotational velocities. The actual accel- 
eration a. l ctrl is a function of the commanded acceleration plus an associated maneuver bias b mvr , scale 
factor s mvr , misalignment /JL mvr , and execution noise T] rnvr terms. 


a ctri ~ T c T (// mvr ) {I + [s mvr \]} T^a ( 


x cmd 


+ * 1 r, 


(27) 


The commanded acceleration is computed using a proportional-derivative (PD) controller 

Kmd = K r { f - r) + K v (v - v) 


(28) 


The reaction control system (RCS) produces the torques required to alter the chaser’s attitude rate. The 
imparted torque is modeled as a function of the commanded torque T c crnd plus a torque bias b rot , scale factor 
s ro t, misalignment fi rot , and noise r\ rot terms. 

T ctrl = T (Vrot) [{! + [Srot\]} Kmd + b roi + V rot } (29) 


The commanded torques are also computed using a PD controller to achieve the desired orientation 

Kmd = Je {Ke (8d c ) + K u (u c c - u> c c )} (30) 

where <5^ c is the 3 -element rotation vector representing the angular difference between the desired attitude 
quaternion q* and the actual attitude q* . The translational and rotational gains in this paper are K r = Kq = 
and K v = = 2£cc n , with natural frequency cc n = 27r/600 rad/s and damping ratio ( = y/2/2. The 

translational and rotational maneuver execution errors are summarized in Tables 4 and 5 respectively. 

DESCENT TO THE ASTEROID SURFACE 

One of the most critical phases of the mission is the descent to the surface of the asteroid. To successfully 
collect the boulder, the capture arms must be centered directly over the object to within 50 cm. The descent 
phase starts 50 m above the surface of the asteroid. The spin rate of the spacecraft is controlled to stay 
above the boulder. Six degree-of-freedom control is applied until the ARV is 20 m above the surface. From 
that point on, control is only applied in the horizontal direction and the gravitational force of the asteroid is 
allowed to pull the ARV to the surface. At 20 m, the downward velocity is required to be 7 cm/s or less. 
Measurements are processed for seven different landmarks which are shown with the trajectory in Figure 
4(a). One landmark is located at the top of the boulder and the other landmarks are 10 and 25 meters away in 
different directions. In this scenario, it takes about 15 minutes to descend from 50 meters above the surface 
to the surface of the asteroid. 
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Bias Initial Uncertainty 

1 

mN 

Scale Factor Initial Uncertainty 

1000 

ppm 

Misalignment Initial Uncertainty 

1 

mrad 

Bias Time Constant 

10000 

sec 

Scale Time Constant 

10000 

sec 

Misalignment Time Constant 

10000 

sec 

Process Noise 

150 

mN/^s 


Table 4. Translational 3 a Maneuver Errors 


Bias Initial Uncertainty 

1 

mN-mm 

Scale Factor Initial Uncertainty 

100 

ppm 

Misalignment Initial Uncertainty 

1 

mrad 

Bias Time Constant 

10000 

sec 

Scale Time Constant 

10000 

sec 

Misalignment Time Constant 

10000 

sec 

Process Noise 

3e-7 

mNm/^/s 


Table 5. Rotational 3a Maneuver Errors 


Control Logic 


The translational and rotational control logic each implement PD controllers as described in Eq. (28) and 
Eq. (30) respectively. The attitude guidance is designed to align the boulder capture mechanism (along the 
ARV axis of symmetry) with the desired line of descent, which is fixed in the asteroid frame. The desired 
angular rate is about the direction of translational acceleration. The translational guidance relies upon a 
pre-planned reference trajectory (relative position and velocity), delivered by the rendezvous and proximity 
operations team at GSFC until the ARV reaches 20 m above the surface. At this point the control logic 
inhibits thruster firings along the direction normal to the asteroid surface to avoid pluming the boulder, letting 
the asteroid’s gravity pull the ARV to the surface. 





(c) Relative 3 a Position Environment Dispersions 


(d) Relative 3a Position Navigation Errors 


Figure 4. Dispersion and Navigation Errors for Descent Phase 
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Results 


A requirement for the boulder capture phase is for the horizontal dispersions to remain within 50 cm. 
Figures 4(b) - 4(d) show the dispersions and navigation errors for this scenario. The dispersions in the altitude 
direction increase drastically after 20 m because no control is applied in that direction but the crosstrack and 
downrange dispersions are well within the requirement, at about 25 cm. 

ENHANCED GRAVITY TRACTOR DEMONSTRATION 

A secondary goal of this mission is to demonstrate planetary defense capabilities. The effect of deflecting 
an asteroid kinematically can be uncertain, especially since many asteroids are rubble piles. The path of 
an asteroid can be deflected by hovering, or in this case maintaining a halo orbit, in front of the asteroid. 
This technique uses the gravitational force of the boulder and the spacecraft to deflect the asteroid. 17 After 
collecting the boulder, the spacecraft will maneuver to a safe distance away from the asteroid and perform 
analyses to better estimate the new mass properties of the spacecraft and boulder. The spacecraft will then 
insert into a halo orbit around the velocity direction of the asteroid, r = 1 km from the center of the asteroid, 
and remain there for six days. The spacecraft will then transition to a halo orbit r = 400 m from the center 
of the asteroid (with a halo orbit radius of r h = 346m) and remain there for 30 days. At this point, the 
Wide FOV Camera and the LiDAR will be blocked by the boulder, so only the Medium and Narrow FOV 
cameras will be available. At this distance the MFOV sensor will be most useful. Because of the location of 
the camera, shown in Figure 2, the camera will only be pointed toward the asteroid during half of the orbit, 
unless the ARV rotates to point the camera in the other direction. The halo orbit is shown in Figure 5(a). 



(a) Gravity Tractor Halo Orbit (b) Dispaced Orbit Gravity Tractor (c) Nominal Trajectory in Inertial Frame 

Figure 5. Gravity Tractor Halo Orbit 


Control Logic 

One option that was investigated for the EGT demonstration was to use only the Solar Electric Propul- 
sion (SEP) System to supply a constant force to counteract the component of gravitational force along the 
asteroid’s velocity vector 17 as shown in Figure 5(b) 


= _ Fo = Gm.im +m,) m (31) 

r z V r / 

where m t ,m c , and are the mass of the asteroid, spacecraft, and boulder respectively. The initial velocity 
for the halo orbit is set so that the centripetal acceleration of the halo orbit equals the vertical component of 
the gravitational force due to the asteroid v z = u, where is the radius of the halo orbit. The radius of 
the halo orbit is also set so that the thrusters do not plume the surface of the asteroid. Another option is to use 
a PD controller to maintain a desired halo orbit. The performance results of both techniques are investigated. 
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(a) Case 1: Dispersed trajectory 




(b) Case 1: Trajectory Dispersions 




(c) Case 1 : Navigation error 



(d) Case 2: Dispersed trajectory (e) Case 2: Trajectory Dispersions 


(f) Case 2: Navigation error 





(g) Case 3: Dispersed trajectory (h) Case 3: Trajectory Dispersions 


(i) Case 3: Navigation error 





(j) Case 4: Dispersed trajectory (k) Case 4: Trajectory Dispersions 


(1) Case 4: Navigation error 


Figure 6. Halo orbit performance results. 


Results 

Several factors such as halo orbit size, orbit maintenance strategy, unmodeled accelerations, measure- 
ment pass duration, frequency of each available measurement pass, camera field-of-view, number of features 
tracked, and resolution influence the overall gravity tractor performance. Although each of these have been 
investigated, this paper will focus on several that influenced the overall vehicle and mission operations, in- 
cluding orbit maintenance strategy, measurement pass duration, and measurement pass frequency. 

The nominal trajectory provided in Figure 5(c) is used to investigate the closed-loop GN&C performance 
with four different test cases: 1) constant thrust profile with continuous measurements, 2) closed-loop PD 
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controller with a 10 minute measurement pass duration once per orbit, 3) closed-loop PD controller with a 
90 minute measurement pass once per orbit, and 4) closed-loop PD controller with a 10 minute measurement 
pass twice per orbit. For these cases, a generic relative sensor is used to highlight important trends, but the 
camera model has also been incorporated to confirm performance levels and show the expected measurement 
pass duration as a function of field-of-view and the number of features. 

Figures 6(a) - 6(c) show the results for Case 1 in which continuous thrust is applied while processing 
continuous measurements. Due to the instability of the gravity tractor orbit, even though the navigation error 
stays tightly bounded and constant thrust is applied, the dispersions grow to almost 2 km. This highlights 
the limitations of using a simple constant thrust halo orbit strategy and emphasizes the need for a feedback 
control implementation for long duration halo orbit maintenance. 

Figures 6(d) - 6(f) show the results for Case 2 which illustrates the advantage of using a closed-loop PD 
controller for maintaining the vehicle on the halo orbit. Although the nominal thrust profiles for the PD 
controller and the constant thrust methodologies are largely the same, the resulting trajectory dispersions are 
drastically reduced from several kilometers to several hundred meters by incorporating a feedback controller. 
For this second scenario, the measurement pass duration was limited to 10 minutes and only one measurement 
pass was allocated per orbit. 

Figures 6(g) - 6(i) show the results for Case 3 where a PD controller is used but the measurement pass 
duration is extended from 10 minutes to 90 minutes. This represents the scenario where either the camera 
FOV is increased or the camera is gimbaled to track landmark features for longer periods. By tracking 
the features for a longer duration, the maximum 3cr RSS trajectory dispersion dropped from 370 m to 70 
m. The improvement comes because the time required to propagate the navigation solution in between 
measurement passes is reduced. Since the trajectory dispersions largely follow the navigation error when 
using a PD controller, the trajectory dispersions are reduced by limiting the error growth due to propagating 
the navigation solution. 

Figures 6(j) - 6(1) show the results for Case 4 which uses a PD controller and only a 10 minute measurement 
pass duration but with the frequency of measurement passes increased from once per orbit to twice per orbital 
period. Although not formally derived at this phase, it is desired that the maximum 3a RSS trajectory 
dispersions remain within 50 m. This accuracy can be accomplished by having short but more frequent 
measurement passes. So even though the measurement duration is only 10 minutes, the maximum 3a RSS 
relative trajectory dispersions are less than 35 m when the frequency is twice per orbit. 

Figure 7 provides a breakdown of the major factors such as initial condition uncertainty, maneuver exe- 
cution errors, process noise (unmodeled accelerations), sensor accuracy, and other environment disturbances 
that contribute to the overall relative position dispersions. This error budget 4 is generated using Case 4 as an 
example because this scenario represents the recommended approach where potentially short measurement 
passes occur multiple times per orbit. In Figure 7, the solid dark (maroon) line represents the 3a RSS relative 
position dispersions given all the potential error sources as previously shown in Figure 6(k). The contribu- 
tion of each error source is then represented with a corresponding line and their RSS value is shown with 
the dashed light (cyan) colored line. These results show the long-term dominating error sources are sensor 
measurement and maneuver execution errors. 


CONCLUSIONS 

These results demonstrate that with the current set of sensors baselined for the ARRM Option B mission 
it is feasible to meet requirements for the boulder capture phase of the mission. Using only a continuous, 
open-loop, thrust profile to maintain the gravity tractor halo orbit, however, is not feasible when accounting 
for unknown perturbations. The instability of the halo orbit requires a closed-loop controller for long term 
orbit maintenance and it is recommended that a second measurement arc be added to decrease the navigation 
error. Future analyses are planned to model thruster locations, throttling, and gimbaling. 
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Figure 7. Relative 3 a Trajectory Dispersion Error Budget (Case 4) 
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